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INTRODUCTION 


The final phase of motion of a satellite equipped with 
a system of gravitational stabilization with a magnetic 
damper Is studied In this work; after completion of the 
transition process It goes Into a stationary state. Then, 
Interaction of dissipative and excitation factors are bal- 
anced In such a way that further evolution does not occur 
and a stationary type of motion Is retained for an Indef- 
initely long time period (for the assumptions made In this 
work). 

A system of gravitational stabilization with a magnetic 
damper Includes a gravitational rod at whose end the magnetic 
damper Is fastened. At a certain stage In motion of the 
satellite, the gravitational rod Is opened to form the re- 
quired ellipsoid of Inertia. 

One of the possible designs for the damper for a gravi- 
tational stabilized satellite is described In reference [1]. 

It Includes a spherical permanent magnet (a damper float) 
placed Inside a spherical cavity In the damper housing. The 
gap between the float and the housing Is filled with a 
viscous fluid In which dissipation of energy occurs during 
relative motion of these bodies. Another dissipation mechan- 
ism is the Foucault currents which occur In the metal housing 
with a relative shift of the float and damper. 

In reference [2] possible periodic motions of this 
system are studied In a plane of a circular polar orbit as- 
suming that the coefficient of damping is fairly small and 
the permanent magnet of tae damper precisely follows the 
vector of directivity of the geomagnetic field. It is pointed 
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out that with these conditions, periodic oscillations and 
periodic rotation with frequency exist commensurate with the 
orbital. In reference [2] asymptotic series which describes 
periodic oscillations produced by stationary solutions of an 
unperturbed system are constructed which converge with fairly 
small values of the coefficient of damping. On the other 
hand, it is clear that when considering the problem posed, 
when the floats of the damper track the direction of the 
force lines of a geomagnetic field, v/lth an Increase in the 
coefficient of damping, the presence of a strong connection 
between the satellite and the damper causes rotation along 
with the magnetic field of the entire satellite. With in- 
finitely large value, the coefficient of damping remains only 
one type of periodic motion, periodic rotation of the satel- 
lite together with the vector of directivity of the magnetic 
field. Consequently, all of the other types of periodic mo- 
tion disappear with certain finite values of the coefficient 
of damping. Computation of this value results actually in 
finding the radius of convergence of the appropriate power 
series . 

In this work, an estimate is made of this value for 
periodic oscillations of a satellite close to axisymmetrlc 
produced by a stable stationary point of an unperturbed 
system and an asymptotic formula is obtained which describes 
this motion up to the maximum value. The asymptotics ob- 
tained of a bifurcated curve (the curve on which the origin 
of pairs of periodic solutions occurs) agrees well with all 
allowable values of the parameters and the results of numer- 
ical computation. 

In all fields of change of the parameters by a numerical 
method, initial data of periodic solutions were obtained 
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both for an oscillation and for a rotating type. For large 
values of the coefficient of damping, where numerical loca- 
tion of initial data of periodic solutions is difficult, 
asymptotic formulas are constructed. 

As is indicated in [2], at certain values introduced into 
the equation of motion of parameters, it has fairly asymptotic- 
ally stable periodic solutions which correspond to different 
stable motions of the satellite with a magnetic damper. Which 
of these motions is realized depends on the initial conditions. 
Therefore, the entire plane of initial values of phase variables 
can break down into those regions so that motion beginning in- 
side one of these regions occurs in a uniform type of periodic 
motion. Each of these regions will be called a field of effect 
of a corresponding periodic motion. 

In this work, the problem of breaking down the space of 
initial data into fields of effect of different types of 
periodic motion will be considered. By numerical integration 
of equations of motion of a satellite for certain values the 
parameters introduced into it, the fields of effect of dif- 
ferent types of periodic motions will be constructed. 
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STUDY OP PERIODIC MOTIONS OP A SATELLITE 
WITH A MAGNETIC DAMPER 

by 

Yu. A. Sadov and A. D. Teterln 


1. Equation of Motion of a Satellite 

In this work, a study Is continued begun In the article 
[2]. As In that article here we propose that a satellite 
moves In a circular polar orbit, the magnetic field of Earth 
Is modeled by a straight central dipole, and the Initial con- 
ditions are selected In such a way that motion occurs In the 
plane of the orbit. Also, It Is assumed that +-he permanent 
magnet of the damper follows precisely the direction of the 
magnetic force line. This assumption Is achieved with ex- 
isting hypotheses as to the values of parameters contained 
in the equation of motion of the float of a magnetic damper, 
namely, when assuming the smallness of the ratio of moments 
of Inertia of the float and the satellite and the smallness 
of the moment of viscous forces acting on the float in com- 
parison with the magnetic moment. 


With these hypotheses, th? equation of motion of a 
satellite relative to the center o' mass will have the fol- 
lowing form 


~ *c<*sin^cos^~S€^- 


\ 0 ^) 

d)) f * 


( 1 ) 


where v - Is the independent variable of width, 

^ - angle of inclination of the largest axis of the 
ellipsoid of Inertia from the local vertical. 


“Numbers in the margin Indicate pagination in the foreign text 


o* - inertia parameter, 

2e - dimensionless coefficient of damping. 




( 2 ) 


Here A, B, C - are the main central moments of inertia of a 
satellite located in increasing magnitude; Kg - is the co- 
efficient of damping; uq - is the orbital angular velocity. 


We note that the right part of equation (1) was apparent ^ 
due to the presence of the damper and introduces excitation 
into the orientation of the satellite. Operation of the sys- 
tem of stabilization can be considered as satisfactory only 
in a case where these perturbations are small. Therefore, 
the dimensionless coefficient of damping e must be fairly 
small. 


For convenience we make a change in equation (1) of the 
variables 

H.l 

( 3 ) 




Then, for describing plane motion of a satellite with a mag- 
netic damper, we find the equation 




(iO 


where differentiation is precisely designated according to T. 
Equation (^) will be the main object of our study. 


As was noted above, parameter e, introduced into equa- 
tion (^), is fairly small. This makes it possible to use, 
in this or another form, a method of a small parameter. For 
this adaptation it is necessary to know the solution of the 
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unperturbed equation which is obtained from the initial when 
e-0. In this case, the unperturbed equation is presented as 
the well known equation of a mathematical pendulum 

• O . (5) 

% 

Equation (5) has the Integral of energy which can be 
written in the form 

• ^ 

^ (ol*^0). ( 6 ) 

Using this integral, it is possible to Integrate equation (5) 
to its end. The general solution of the equation is expressed 
by elliptical Jacobi functions [3]. 


Where 0<v<l, equation (5) describes oscillation motion 
with the period 




A'* V , 


(7) 




where K(k) - is the full elliptical Integral of the first order. 


When v>l, equation (5) has a periodic solution which de- 
scribes rotation with the period 


55 . , 



( 8 ) 


The case v*l corresponds to motion along the separatrlx 
separated from the phase plane ( 0 , 0 ), the field of oscilla- 
tions from the field of rotation. 


Besides these solutions, equation (5) has two stationary 
solutions 

QmO, 9*0 (V*0) , 

e>jT, 9‘0 (v-j) . 

3 




Equalities (10) define the trough special point of unperturbed 
motion. Tlie second special point is the center which is de- 
scribed by equation (9)* 

Formulas (6)-(8) lose their meaning when a«0 (an axlsym- 
metrlc satellite). In this casei equation (5) describes a 
uniform rotation 

The picture of the trajectory, as is seen from equation 
(5), is periodic according to e, and therefore the phase plane 
(0, 0) can reduce to a cylinder (-ir<0£ir; -•<e<+»). All of the 
phase trajectories on this cylinder except for the separatrlx 
are closed and periodic according to T, 

2. Classes of Solutions Generated 

Periodic solutions are of particular Interest for studying 
equation (^1). Equation (^) is periodic according to the vari- 
ables T and 0, that is, its form does not change with substi- 
tutions 

T-^T+Pir , 

0-*-0+2lt. 

Therefore, periodic solutions of this equation must satisfy 
the condition 

9(Z^. jrm)-9(T)*2jrp (,n- ... . ^ ( 12 ) 

The trajectory which satisfies this condition for m 
periods according to T circles the phase cylinder (0, 0) p 
times and is closed going to the initial point. When e-^0, 
such a trajectory changes to a certain integral curve of un- 
perturbed motion which is called generated for the solution 
considered. Due to the continuous dependence of the solution 
on the parameter c, the generated solution also satisfies 


condition (12). This makes It possible to separate^ among 
the solutions of an unperturbed system, those which can be 
generated for periodic solutions of a perturbed system. 


When a«0 from equations (11) and (12) we find that on 
the generated trajectories one must fulfill the condition 





(13) 


If aWO, then among the generated solutions one finds station- 
ary solutions (9)» (10) and also oscillatory (7) and rota- 
tional (8). 


Let us consider oscillation solutions which are defined 
by formula (7). For these, the relationship following is 
fulfilled 

0(t+Ti) ■ 0(t) 


and, consequently, 

e(T*nTf ) • 9(T) , > 7 - 3 , . . . 


(liO 


Comparing (l^l) and (12) we find that on oscillation trajec 
tories, (12) is fulfilled when p“0. If only 





(15) 


For rotational solutions from formula (8) we have 


S(z*r,).e(T)tsx , 

or 

6(Z*nr^)m ... 

In this way, condition (12) is fulfilled, if 


( 16 ) 


m 


(17) 


Then, in formula (12) p*in, where the sign for n is selected 
in accordance with the direction of rotation of unperturbed 
motion. 

In reference [2] it is pointed out that periodic solu- 
tions do not exist for equation (4), generated by oscilla- 
tion solutions (7) of an unperturbed system. Also periodic 
solutions do not exist generated by rotations (8) toward the 
side opposite the direction of rotation of the magnetic field. 
Rotation (8) in a straight direction can generate slow (with 
a frequency no larger than the orbital, m^^n) nondamped rota- 
tions of the satellite. With this rotation at an orbital 
frequency (m«n) it is possible with all physically allowable /13 
values of the parameters, (0<o2^3» 0_^c<*»). Periodic rota- 
tions of other types have a very narrow field of existence 
in the space of parameters of o, e and, accordingly, a low 
probability of falling into this condition. Stationary so- 
lutions of an unperturbed system are generated with stable 
(9) and unstable (10) oscillations. When a»0 (an axisymmetrlc 
satellite) for any e>0, a single parameter family of stable 
periodic trajectories exists, corresponding to m«p*l: 

e.T.sf: 

TTt Hi** £** H* 

3. Construction of Periodic Solutions of the Equation of 
Motion of a Satellite Close to Axisyminetric . Generated 
by Stationary Solutions of an Unperturbed Equation and 
an Evaluation of the Field of their Existence in the 
Space of Parameters (c. <t ) . 

Having u.'3ed thje smallest values of parameter c, one can 
construct the asymptotic formulas wli’ jh describe periodic 
motion of a satellite, generated by stable (9) and unstable 
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(10) stationary points of unperturbed equation (5). As Is 
pointed out in reference [2], they have the following form 



(18) 

(19) 


for all physically allowable values of the Inertia parameter /l^ 
a(0<a^/3) • 

With the proposed posing of the problem (when the float 
of the damper follows the direction of the force lines of 
the magnetic field) 2ir - periodic oscillations (l8) and ( 19 ) 
will exist not for all values of a dimensionless coefficlenl 
of damping c, but only for 0<c<c*, where the value of de- 
pends on a. Actually, for fairly large values of c, the 
presence of a strong connection between the satellite and the 
damper float causes rotation along with the magnetic field of 
the whole satellite. Rotation of m:n*l:l Is a unique type of 
periodic motion In this case. 


Numerical values of e,, result actually In finding the 
radius of convergence of series (l8) and (19)* This problem, 
as a rule. Is extremely cumbersome, that is. It requires com- 
putation of the common member of the series (l8) and (19) 
which In the case of nonlinear equations presents considerable 
difficulty. Therefore, we will attempt to attain certain 
estimates of the value of eg. For this, actually we will con- 
struct a periodic solution generated by the stable position 
of equilibrium (9) In a case when the Inercia parameter a Is 
small. Then one can expect that the critical value of e* also will 
be small and one' can use asymptotic methods right up to this 
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value. Here, however, it is impossible to use asymptotics 
of (l8), because the amplitude of oscillation when c is 
close to is large. On the other hand, angular velocity 
0, as one sees from the integral of energy in unperturbed 
motion (6) remains a small value on the order of o. 


Therefore we assume Q^at-p 


Then 




( 20 ) 


( 21 ) 


/15 


We assume 




( 22 ) 


and rewrite system (20), (21) In the form 

6 » Hp % ( 


( 23 ) 


Finding the solution of this sytem in the form of the series 
according to the power fz 

pmp* Vepf * spt * *■ 


' I • • • 


( 24 ) 


we find the following equations for sequential approximations 

' 0^‘O, Po'O-. 

k-xp., ( 25 , 

dt COS t,-pj i 

p^M--j9gCos0c^* -^QfSinQo-pf ; • • • 

Using the solution of the equations of a zero approximation 

0 ^C,-const , • ( 26 ) 
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we find in the first approximation 

Sis 


where 


9t^Md,T*C, , - ffSMC/Z^: 


(27) 


/16 


( 28 ) 


For this, so that solution (27) would correspond to periodic 
oscillation, it is necessary and adequate that 


Dg* o 


9 



(29) 


The latter condition can be fulfilled only when k^2. 

For the calculation of (22) this result can be formulated as: 
when 

(30) 

does not exist of periodic solutions of equation (^), gener- 
ated by a stable stationary point. The value of e*, defined 
by formula (30) gives an estimate above the boundary of 
of the field of existence of periodic solutions of the type 
considered (e*<e'). 


The constants Cj and entering into (27) are determined 
from the condition of periodicity of the solution of equations 
of the second approximation 


0 ^. 3 . 

p^‘--^C,cosCo*Di, . 


(31) 


From this we find 




f 



(32) 
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In a completely analogous way, the conditions of period- /17 
icity for all sequential approximations give an equation for 
determining the additive constants of the preceding approxi- 
mation. These equations are linear and solved when cosGof^O; 
therefore, construction of the series (2^) is possible for 
all e, Which satisfy the inequality 


The dependence of g' on a is graphically shown in Figure 1. 

^ • Here, along the axis of 

t theaJs's^slssa, the param- 



Flgure 1 


eter a=a2/4 Is applied 
and e’(a) is expressed 
by a dashed straight 
line. The solid line in- 
dicates the boundary of 
the field of existence 
of the solutions con- 
sidered, obtained using 
numerical calculations 
described in the follow- 
ing paragraph. 

A total of the con- 
stants C 2 and D 2 in form- 
ula ( 31 ) gives us 




s -0.73ZS 


(3iJ) 


Consequently, with small values of a, the periodic solutions 
sought can be written in the form of a series 






(35) 
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^ . Numerical Finding of Periodic Solutions 

In a broad range of values of the parameters a and e, 
periodic solutions are found by numerical Integration of equa- 
tion (^). The problem amounted to a selection of the Initial 
conditions 0, 0 when T=0 so that at the moment T»2ir the fol- 
lowing relationship is fulfilled 

^Q(o) ^ Q(^jr) • 0(0 ) . ( 36 ) ^ 

In this way, periodic oscillations were calculated, generated 
by stable and unstable stationary points and also periodic ro- 
tations of the type m:'’=l:l. Computation encompassed a range 
of values from 0 to 0.750. The latter number corresponds 

to the boundary of physically allowable values of the Inertia 
parameter. 

The values of the initial conditions of integration ob- 
tained from computation 0^°^=0(o) and 0(®^=0(o) depend on the 
two parameters a and e 

and are shown in Figure 2. In this drawing, three series of 
curves are visible. Each of the curves is the geometrical lo- 
cation of the initial conditions of periodic solutions v;lth a 
fixed value of the parameter a and a change in e from 0 to 
In each series, it is expressed on the three curves correspond- 
ing to a=0.3; 0.6 and 0.75. 

Curves in the upper part of the drawing when 0(°^>1 cor- 
respond to rotation, then when 0^°^>O, stable rotations are 
obtained; when 0(°)<O - they are unstable. Such motions ex- 
ist with all values of e and therefore the curves of each 
series begins when e=0 and are completed in general for each 
of the series at a point corresponding to £-♦■». In this maxi- 
mum case, the satellite like the float of the damper rotates 
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/31 



i 0(0) 


Figure 2 


along with the force lines of the magnetic field and such mo- 
tion does not depend on gravitational moments. The numerical /19 
location of the solutions with large values of e are diffi- 
cult but in this case one can construct asymptotic formulas: 





- for unstable 
rotation. 

(38) 


- for stable 

(39) 


rotation 


formulas p is constant 



p- j{3tn3-t3en2)^ S,33B $ . 


(40) 


REPHODUdBlLirV OF 
page is poor 


For obtaining asymptotic formulas (38), (39) we rewrite 
equation (^) in the form 




S~ScasT £ » 


(41) 


Solution of equation (4l) will be found in the form of a 

series: _ 

/ 


rv ... 


(42) 


By substituting (42) in equation (4l) and equalizing the co- 
efficients of uniform powers of the parameter e, we fltid that 
the function 0j^(T), 1=0, 1,2,... must satisfy the following 
equations: 




•“ S~ 3COST 

G,=-9^-asinQ„ , 


(43) 

(44) 


Q^^~df‘Cf6tcos9o » 

s-d^-a(9j,cosdo-^ sindo) , 


(45) 

(46) 


/29 


Then for any value of T one must satis’fy the condition 


tlon 


&^(r*axn)~9oCt)*SJrn, 9i(T*SJrff)»9^-W, 

Integrating equation (43), we have in a zero approxim^- 

. (i^8) 


The zero approximation describes the motion of the satellite 
when e=“, that is, when it rotates along with the force lines 
of the magnetic field. 


Prom equation (44) 


we find the first approximation 
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t 


. iaasC. &, , 

S-3C0ST 3 ^ & 

*c,sinC^fr.lAi,itg(etf^pC, . 


( 19 ) 


From the condition of periodicity (47) it follows that 


that is. 


sinC^»0, 




(50) 


The first value of the constant corresponds to unstable 
rotation, and the second to stable. 


The second approximation we find by integrating equation /21 
(^5). sirtT 


Q t.gacosC. ^ 

(S~3costf ^S‘3cosT 

.oCa>sC,[fT-fA>cig(ei^§)] • 

3 JS-3COSX S-3co$x * 


(51) 


The condition of periodicity (4?) for 02 (T) takes on the fol- 
lowing form 


23t 

3 ^ i S'Seosx S‘3cosx 

0 


(52) 


Having used formulas 4.413 (3-4) [5], we find 

Cf’CrpcosCo , 

v/here p is determined by formula (40). 


(53) 
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The constant C 2 In (51) we will define from the condition 
of periodicity for the third approximation. It appears equal 
to zero: 


» O . 


(54) 


Thus, for adequately large values of e, equation (4) has 
two 2n-periodic solutions of a roatating type (unstable and 
stable, respectively) which according to formulas (42)-(54) 
can be presented in the form 




f*/ (S’ 


/SsiffT 


I S-ScosT 
'Scasry S’ScosT 




±r^ 

3 


Scosx-S ^ S’Seosx 


B(T)* 


3 { S-Seosx 

O 


& 

<fSt/rZ 

S-Scosz £ S-Seosz 


dx 


(s-seosT'^’)* 


, ± f f2(6’Seo sZ-3eos^) Sa(ScosT’S) 
(S-Scoszy * (S’3cffszy 


-a 


tSeasz-3 

S-Scffsz 




f_ J fSsinZ SasinZ 
TUCS’Seoszy' S-3COST ^ 


Q(Xh 


, ^Fseosx-3 S’3eosx ^j\) n/f\ 
3 1 T’ Scosx 

S-3eosX € S-3COSZ [s-3eosz J 
± riZ(6-Scos7-3CffS*T) 8a(ScosZ~3) 
€*L (S’Seoszy 


-a 


« ScosT-3 


■fp-r^ 


(S’Scoszy 


I 


S-SeosT 


When T-0, from (55) and (56) we find asymptotic formulas (38) 
and (39). 

In the lower part of Figure 2 when there is a /23 

series of curves corresponding to oscillation motion of the 
satellite. Each curve consists of two branches — stable and 
unstable. These branches occur when e»0 respectively from 
the stable and unstable stationary points of an unperturbed 
system. With an Increase in e, the points shown for stable 
and unstable periodic solutions move along their branches 
and when e»e^(a) they merge at the branching point. With large 
values of e, 2iT-perlodlc oscillation solutions do not exist. 

The bifurcated curve on which merging of two periodic solu- 
tions occurs is shown in Figure 2 by the dashed line. The 
critical values of ej, depending on a are shown by the solid 
line in Figure 1. This curve, for all physically possible 
values of the parameter a coincide fairly well with straight 
line 

e-a, (57) 

defined by asymptotics (27). 

The fields of existence of two types of periodic solu- 
tions according to the parameter e are considerably narrower. 

For example, when (a0.6l6) ^Tr-periodic rotations cor- 
responding to m;n«2:l, it can be successfully obtained numer- 
ically only when e<0.0035. 

5 . Fields of Effect of Different Types of Periodic Motion 

As one sees from reference [2], besides the periodic so- 
lutions studied above, with certain values of the parameters 
a and E, equation (^) has other types of asymptotically 
stable periodic solutions which correspond to different 
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established motion of a satellite with a magnetic damper. /2k 

Which of these motions is realized depends on the initial con- 
ditions, Therefore, the entire plane of the initial values of 
Q can be divided In these fields so that motion, beginning 
inside one of the fields results in a uniform type of periodic 
motion. 


Let us determine the transformation of the phase plane of 
6, S for a period like this transformation. In which a given 
point 0, 0 Is the point (0(2u), 0(2ir)) of the trajectory of 
equation (h) which occurs at the moment T*0 through (0(0), 

0(0)). By virtue of the 2ir-periodiclty of equation {k) ac- 
cording to the T properties of the trajectory, corresponding 
to point (0(0), 0(0)), like the initial, are invariant with 
such a transformation. The invariant. In particular, is broken 
down in the field of effect. The 27r-periodic solution of 
equation (4) corresponds to immovable points of the transfor- 
mation considered. It is clear that asymptotically stable 
stationary points belong to the field of effect on the periodic 
solution determined by them. Invariant curves which pass through 
unstable stationary points are the separatrixes of a given 
transformation; they divide the field with qualitatively dif- 
ferent behavior of the invariant curves. In particular, the 
sum of the separatrixes divide the field of effect. Therefore, 
the problem of finding the fields of effect leads to calcula- 
tion of the separatrlx of transformation of the phase space 
caused by a shift in the period along the trajectory of 
differential equation (4). 

As one sees from what has been presented above, such a 
transformation for any fixed values of the parameters a=a^/4 
and e (0<ai0.75; 0<e<cu(a)) has four stationary points (Fig- 
ure 2): two stable and two unstable. Computation of the 




separatrlx began In the environment of unstable stationary /25 
points* Close to such a point, the transformation of the phase 
plane Is close to linear. The Intrinsic vectors of this linear 
transformation Indicate the direction of the branches of the 
separatrlx, which pass through a given stationary point. The 
Intrinsic vector which has an Intrinsic value In absolute size 
larger than one corresponds to outgoing branches of the separ- 
atrlx and the Intrinsic vector with Intrinsic value according 
to the model smaller than one corresponds to the Incoming 
branch. For constructing the separatrlx, the Initial condi- 
tions are selected lying In the environment of unstable sta- 
tionary points In a direction determined by the Intrinsic 
vector and numerical integration is made of the system with 
these initial conditions. The points of the trajectory cor- 
responding to the moment of time T»2kir, ...» belongs 

to the appropriate separatrlx. For this, In order to obtain 
in this way the incoming branches, integration must be done 
in a direction of decrease of the independent variable (T«2kt, 
k“l,-2,...). Obtaining the sequence of points when moving away 
from the Initial stationary point becomes very rare and does 
not give the correct concept of the behavior of the separatrlx. 

For obtaining a more detailed picture, one must follow the In- 
dicated procedure more than once, selecting as the Initial dif- 
ference of the point one or another branch In the environment 
of the unstable stationary point. 

Figure 3 shows the breakdown of the phase plane obtained 
in this way corresponding to a*0.750, c«0.1. In this case, 
as has already been Indicated, we have four periodic solutions: 
two stable and two unstable. In Figure 3» the four stationary 
points correspond to it. At each of the unstable stationary 
points, the separatrlxes converge along four branches: two 

Incoming and two outgoing. These curves break down the phase 
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cylinder In four parts. One of these parts shaded In Figure 
3 forr,^ the field of effect of the rotating periodic motion, ^6 
the t..i’ee others belong to the field of effect of oscillation 
motion. One of the latter three sections corresponds to the 
trajectory beginning in the field of reverse motion with nega- 
tive values of S. The trajectories from the remaining two 
sections begin in the field of forward rotation and end with 
periodic oscillations; they differ from each other in that 
the trajectory from one section (in Figure 3 it corresponds 
to the broader nonshaded band in the upper part of the draw- 
ing) separate only rotation around the phase cylinder in com- 
parison with the other ti ajectories . 



Figure 3. (a»0.75; e»0.1) 

From Figure 3 it is apparent that the field of effect 
of a rotational periodic condition Is fairly narrow and Is 
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located as a whole in the field of forward rotation when S>1. 
Therefore, one of the methods to avoid capture of the satel- 
lite in rotation consists of imparting to it a certain initial 
twist in rotation opposite the direction of orbital motion. 

When increasing the coefficient of damping c, the field 
of effect of the rotational condition is expanded and at cer- 
tain values of c the separatrlx going away from the point 
corresponding to unstable oscillation falls at a point cor- 
responding to unstable rotation. After this, the picture of 
breakdown of phase space in the field of effect changes 
qualitatively. In the field of effect of the stationary ro- 
tational condition there is a band passing into the field of 
reverse rotation. Thus, for example, when a«it/2 and e*0.^, 
one can find the initial conditions in the field of reverse 
rotation leading to stationary rotation of the satellite. A 
further Increase in c results in a larger contraction of the 
field of effect of the oscillation condition which disappears 
when c"£j|(a) and rotation along with the magnetic field re- 
mains the single type of established motion. 

Conclusions /27 

In this work, periodic motion of a gravitationally 
stabilized satellite with a magnetic- dinner In the plane of 
a circular polar orbit is studied. Inthesystem considered, 
besides the small established rotations of the satellite 
relative to the oriented position, also small (•..'1th frequency, 
not large orbital) undampcc: rotations are also possible. 

Then, rotation with the orbital frequency is possible for all 
physical values of the inertia parameter and for all values 
of the coefficient of damping. The periodic rotations of 
other types have a very narrow field of effect and, correspond- 
ingly, a low probability of falling into this condition. 
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With large values of the coefficient of damping due to 
the strong connection between the satellite and the float 
of the magnetic damper, a constant magnet which precisely 
tracks the force lines of the magnetic field, the satellite 
also oeglns to rotate along with the float. Therefore, per- 
iodic oscillations of the satellite around the oriented posi- 
tion exist only for values of the damping coefficient which 
do not exceed a certain maximum value depending on the value 
of the inertia parameter. 

In the work, an evaluation is obtained higher than the 
maximum values of the coefficient of damping for a satellite 
close to axlsymmetric (at low values of the inertia parameter). 

The results of numerical calculation showed that the value 
obtained corresponds well to the actuality for all physically 
possible values of the inertia parameter. Prom the evalua- 
tion it follows that for coefficients of damping .:>0.75, 
periodic oscillations of the satellite do not exlc'* nor at 
such physical values of the inertia parameter as a(0<a;^0. 75) . 

At certain values of the parameters a and c, the equa- /2S 
tlon of motion of a satellite has several a^^ymptotically 
stable periodic sol'itlons which correspond to different 
stablllz«?d motions of tht- satellite with a magnetic Oamper. 

Which of these motions Is realized depends on the Initial 
conditions. Therefore, the entire plane of initial values 
of phase variables can break down in the field of effect 
of different types of periodic motloris, that is. Into those 
fields where notion beginning Inside one or anothei' field re- 
sults in a uniform type of periodic motion. In the work, a 
breakdown of the phase plane In fields of effect for values 
of the parameter a*O.790f c«0.1 Is obtained. In this case, 
the equation of motion has two asymptotically stable 2it- 
perlodlc solutions; one is oscillatory and the other Is 
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rotational. As Is seen from the breakdown, the field of 
effect of the rotational periodic condition is fairly narrow 
and is located as a whole within the field of forward rota- 
tions. Therefore, one of the methods to avoid capture of the 
satellite in rotation is to impart to it a certain initial 
twist in a direction opposite the direction of orbital motion. 
With an increase in the coefficient of damping e, the field 
of effect of the rotational periodic conditions is increased 
and the oscillatory is decreased and when e = e,j(a), the field 
of effect of the oscillation periodic condition disappears. 
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